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We report on recent work on the determination of low-energy constants describing AS = 1 weak transitions, 
in order to investigate the origins of the AI = 1/2 rule. We focus on numerical techniques designed to enhance 
the statistical signal in three-point correlation functions computed with overlap fermions near the chiral limit. 



1. INTRODUCTION 

Much recent activity in lattice QCD has 
focused on the determination of low-energy 
constants (LECs), which parameterise non- 
perturbative physics in an effective description at 
low energies. A promising approach involves com- 
puting in the e-regime, the kinematical region of 
arbitrarily small quark masses m in a finite vol- 
ume V such that rnEU < 1. Here discretisations 
that preserve chiral symmetry at non-zero lattice 
spacing should be used so that small masses and 
momenta can be reached safely. However, sim- 
ulations in the e-regime are expensive: owing to 
spontaneous chiral symmetry breaking the spec- 
trum of the Dirac operator becomes arbitrarily 
dense near the origin, and thus the operator is 
ill-conditioned. Furthermore, one observes large 
statistical fluctuations ("spikes") in correlation 
functions of local currents |1I2| for small masses. 
A number of efficient numerical tools have been 
developed |3I4I2| . in which the low modes of the 
Dirac operator are treated exactly. In ref. [2j we 
proposed a method, dubbed "low-mode averag- 
ing" (LMA), which was shown to eliminate the 
spikes in two-point correlation functions for quark 
masses m < 1/SV. Here we report on the exten- 
sion of LMA to three-point correlation functions, 



in order to study non-leptonic kaon decays. 

2. LMA REVISITED 

We consider the massive Neuberger operator jS] 

D m = (1 - \am)D + m, a = a/(l + s), (1) 

where \s\ < 1 is a free parameter and the massless 
Neuberger operator, D, satisfies the Ginsparg- 
Wilson relation. We are interested in correlation 
functions of the left-handed flavour current 

W«? = (U^), P± = |(l±75), (2) 

where ip = (1 — aD/2)ijj, and a, (3 are flavour in- 
dices. Using the left-handed current has the ad- 
vantage that zero modes of the Dirac operator are 
projected out in correlation functions. However, 
modes corresponding to small non-zero eigenval- 
ues can still cause large statistical fluctuations. 
Let S(x, y) denote the quark propagator from y 
to x. LMA works by separating m ow low-lying 
modes and using the spectral representation of 
the propagator. In this way the quark propaga- 
tor is written as 
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S(x,y) 



efc(z) g e k (y) 
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S h (x,y), 



(3) 



*Talk presented by H. Wittig at Lattice 2004, Fermilab, 
21-26 June 2004. 



where S h is the propagator in the orthogonal 
complement of the subspace spanned by the n\ ow 
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lowest modes. The vector e k is given by 

e fe = P a u k + P^ IJ DP IJ u k , (4) 

where — <r denotes the chirality of the sector con- 
taining zero modes (if any), and u k is an approx- 
imate eigenmode of D^Dm. 

P^D^DmPo-Uk = a k u k + rfc, (u fc , n) = 0. (5) 

After inserting the rhs. of eq. J3J into the expres- 
sion for the correlator C(t) = ^2g{Jo(x)Jo(0)) , 
one picks up three contributions: 

C(t) = C u (t) + C h \t) + C hh {t). (6) 

As explained in the statistical signal for the 
correlator is enhanced by exploiting translational 
invariance in C 1 (t) and C hl (t) , so that these con- 
tributions are sampled over many different source 
points. Applications of LMA to the calculation of 
in the e- regime can be found in 0. 

3. K — ► 7r7r WITH ACTIVE CHARM 

In a recent paper [H] we proposed a strategy to 
investigate the origins of the AI — 1/2 rule in 
K — > 7T7T decays, i.e. the observed enhancement 
of the transition amplitude for a two-pion final 
state with isospin 0, | ^4q | / 1 ^2 1 = 22.1. 

The goal of our programme is to quantify sep- 
arately the contributions to the AI =1/2 rule 
from physics at the charm quark mass scale and 
from "intrinsic" QCD effects due to soft-gluon ex- 
change. A crucial part of our strategy is to keep 
an active charm quark, so that the theory has an 
approximate SU(4)l x SU(4)r chiral symmetry. 
The AS = 1 effective weak Hamiltonian after the 
Operator Product Expansion contains the four- 
quark operators Qf and Qj~ which are given by 

Qf = (s^P-u){u^P-d)±(s^P-d)(u^P^u) 
-(u-c). (7) 

Qf and Q± are singlets under SU(4)r and trans- 
form under irreducible representations of SU(4)l 
of dimensions 84 and 20, respectively. 

In Chiral Perturbation Theory (ChPT), the ra- 
tio of amplitudes jvlol/l^l is, at leading order, 
related to a ratio of LECs via 

14)1 _ i n 3<?n m 



Here, the LECs gf multiply the counterparts of 
the four-quark operators Qf in the low-energy 
theory. They can be determined by studying 
suitable correlation functions of Qf in QCD and 
matching them to the corresponding expressions 
in ChPT. The fact that in the e-regime this 
matching can be performed at NLO without the 
presence of additional interaction terms with un- 
known coefficients 6 , makes this an attractive 
setting for carrying out our programme. 

The intrinsic QCD contributions to the en- 
hancement can be isolated by determining <?j~ / gf 
in the theory with m u = rrid = m s = m c . Then, 
by studying the ratio of LECs as m c departs from 
the mass-degenerate limit, the specific contribu- 
tion of the charm quark to the AI = 1/2 rule can 
be investigated in detail. 

The complicated renormalisation and mixing 
patterns of four-fermion operators usually en- 
countered in lattice formulations can be avoided 
through the use of Ginsparg- Wilson fcrmions. In- 
deed, the use of the modified fields tp = (1 — 
aD/2)ip in eq. guarantees that the renormali- 
sation and mixing of Qf are like in the continuum 
theory [Sj. In particular, no mixings with lower- 
dimensional operators can occur. For full details 
we refer the reader to ref. jo]. An investigation 
of the effects of the charm quark based on ChPT 
has been published in 

We now report on our attempts to extract the 
ratio gi/gf in the SU(4)-symmetric theory. In 
this case only "figure-8" -graphs must be consid- 
ered, since "eye" -graphs cancel exactly for m c = 
m u . We define the following three-point functions 
of Qf and the left-handed flavour current: 

c£(xo,vo) = E ( x )Uu [Qf (Q)] [Jo (y)]us ) ■ (9) 

The relation between g^jgf and the ratio of cor- 
relators is given by 

— H{x ,y ) = - T — r , (10) 

where the factor H(xo,yo) has been computed in 
ChPT in the e-regime at NLO 0. The factors 
kf are Wilson coefficients, and it is assumed that 
C-f jCf is rcnormalised in some scheme. From 
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Figure 1. The fitted value of the ratio Cf/Cf 
plotted versus mSV for n\ ovl = 5. We have set 
£ = (250 MeV) 3 . 

eq. Q10[) one sees that the ratio of correlators is 
directly proportional to the ratio of LECs. 

Our main results were obtained from 638 con- 
figurations on a lattice of size 12 3 ■ 30 at j3 = 
5.8485, using the 5 lowest modes in the imple- 
mentation of LMA for figure-8 graphs. In Fig.^ 
we plot the ratio C-f /C^ as a function of the 
quark mass in units of EV. It is seen that Cf/Cf 
increases as the quark mass is tuned towards 
smaller values, which, at first sight, might be in- 
terpreted as an enhancement in the ratio of LECs 
from QCD effects alone. However, the quark 
masses used in the simulations are too heavy to 
lie inside the e-regime, so that the computed NLO 
matching factor H(xo,yo) cannot be applied to 
yield g{ j g\ . For smaller masses the statistical 
signal for the ratio of three-point functions was 
lost even with LMA. 

In order to demonstrate the feasibility of our 
strategy one has to find ways of obtaining a strong 
signal in the e-regime. We have thus investigated 
the quality of the signal for three-point functions 
using LMA with a larger number of low modes 
[Hj . Preliminary results on a lattice of size 8 3 • 20 
at j3 = 5.8 are shown in Fig. [3 where we have 
plotted the bootstrap distributions for the ratio 
of three-point functions C-f /Cf. The plot shows 
that using m ow = 8 hardly improves the signal at 
all, while for n\ ow = 20 the distribution is much 
narrower, so that the statistical signal is indeed 
enhanced. We conclude that LMA leads to an 
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Figure 2. Bootstrap distributions obtained for 
mHV = 0.917 at [3 = 5.8 on a 8 3 • 20 lattice. 

improvement also for three-point functions, but 
only if niow is increased relative to the previously 
studied case of two-point correlators. Obviously, 
using larger values of n\ ow leads to a computa- 
tional overhead, which, however, is outweighed 
by far through the observed enhancement of the 
signal. We are currently running on larger spatial 
volumes in order to corroborate these findings. 
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